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Abstract. In this paper, we describe an algorithm for computing al- 
gebraic modular forms on compact inner forms of GSp4 over totally 
real number fields. By analogues of the Jacquet-Langlands correspon- 
dence for GL2, this algorithm in fact computes Hecke eigensystems of 
Hilbert-Siegel modular forms of genus 2. We give some examples of such 
eigensystems over ' 



Introduction 

Let F be a totally real number field and D/F a quaternion algebra over 
F that is ramified at all infinity places of F. Let G/Q be the algebraic group 
whose Q-rational points are given by the (quaternionic) unitary similitude 
group GU2(-D). Then, G is an inner form of G' := Resj7'/(Q(GSp4) such 
that G(M) is compact modulo its center. The Langlands philosohpy pre- 
dicts that there is a transfer map between automorphic forms (or algebraic 
modular forms in the sense of Gross [11]) on G and Hilbert-Siegel modu- 
lar forms on G'/Q. In fact, the existence of such a transfer map was first 
conjectured by Ihara, Hashimoto-Ibukiyama, and Ibukiyama [16 ^ 113 ^ [T5] for 

= Q, as an analogue in genus 2 of the Jacquet-Langlands correspondence 
for GL2 or, more precisely, the Eichler theta correspondence. And recently, 
Sorensen [23] established that such a correspondence does indeed exist when 
the degree of F is even. 

Under the assumption of the existence of a Jacquet-Langlands correspon- 
dence between G and G' , Cunningham-Dembele [5j give an algorithm for 
the computation of Hecke eigensystems of Hilbert-Siegel modular forms of 
genus 2 over Q(\/5). Here, we extend the algorithm in that paper to any ar- 
bitrary totally real number field F. Our approach borrows from Hashimoto- 
Ibukiyama and Ibukiyama [13^ [T5] who provided numerical evidences in sup- 
port of their conjecture by explicit computations with Brandt matrices. And 
also from Lansky-Pollack [17] who computed examples of algebraic modular 
forms on G for F = Q and determined their Satake parameters. (We note, 
in passing, that the computations in [TTj were used to predict the existence 
of symmetric cube liftings from PGL2 to PGSp4, a result which has now 
been proved by Ramakrishnan-Shahidi [18j.) 
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By work of Ghitza [10], all Hecke eigenvalues systems of modp Siegel 
modular forms come from Hecke eigenvalues systems computed from the 
compact form G obtained from the quaternion algebra D/Q ramified at oo 
and p > 3 only. His result combined with our algorithm provide a mean to 
gather numerical evidence on the "weight part" of the analogue of the Serre 
conjecture for GSp4 (Herzig and Tilouine [Hj), and we hope to do that. 
With some changes, this algorithm should also works for other algebraic 
groups which satisfy the conditions in Gross [111 Proposition 1.14], especially 
compact forms of U(3). And we also hope to use it for this latter group 
in order to investigate a conjecture of Bellaiche-Graftieaux [Ij which is an 
analogue of the so-called Ihara Lemma. 

The outline of the paper is as follows. In Section [H we give a brief review 
of compact inner forms G of G' and the theory of lattices on integral models 
of G. In Section [21 we recall the theory of algebraic modular forms on G. 
Then in Section [3l we give some details regarding the implementation of the 
algorithm. Finally, in Section U] we give numerical examples over the real 
quadratic field Q(\/2)- 
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1. Inner forms and integral models of GSp4 

Throughout this paper, we let be a totally real number field of degree g 
and Op its ring of integers. For every real embedding v : F ^M. and every 
element a £ F, we let a„ := v{a) be the image of a under v. (We recall that 
an element o G -F is said to be totally positive if a^, > for all v \ oo.) Let D 
be a quaternion algebra over F which is ramified at all infinite places; we let 
S be the set of finite primes at which D is ramified. We denote by i— > n and 
by nr : — > F the involution and the reduced norm of D, respectively. We 
fix a maximal order Oo in F>. We choose a finite extension E/F, contained 
in C, and such that there is a splitting isomorphism j : D 0q E ~ M2(-E')^. 
For any prime p C Op, we denote the completions of F and Oi;' at p by Fp 
and Opp, respectively. Similarly, we denote the completions of D and Od 
at p by Dp and Od^- For any finite prime p ^ S, we fix an isomorphism 
~ M2(C'fJ and extend it to Dp ~ M2(Fp). 

1.1. A review on O^j-lattices in D^. Here, we recall some results about 
lattices in and supplement them. We consider with its left -D-module 
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structure. The endomorphism ring of L>2 is M2{D), the ring of 2 by 2 
matrices with entries in D. The involution of D induced an involution on 
M2{D) by action on the entries. Let A be a matrix in M2{D). We define 
the determinant of A as follows. We choose an extension E/F which splits 
D and let det(^) = det£;(A (g) 1) in M2{D Of -E) ~ M4(^). It is not hard 
to see that det(^) does not depend on the choice of splitting field. We say 
that A is hermitian ii A = A^. In that case, 

A = with s,t e F, r e D, 

and we define detniA) := st — nr(r). Then, we see that det(^) = det/)(A)^. 
A quaternion hermitian form on D is a map H : D^x ^ D such tat 

H{u, v) = H{v, u) and H{au, hv) = aH{u, v)b y a,b D, u,v D^. 

By choosing a basis, we can write H as 

H{u,v) := uAv\ u,veD'^, 

where ^ is a hermitian matrix. We say that A is totally positive and write 
A > a the form (u, v) i— > uAv^ is totally positive definite. 

Prom now on, we fix such a totally positive form H on and let A be 
the associated matrix. A left O/j-lattice L C is an Oi^-lattice in 
which is also a left O/j-module. As in Shimura [22, Section 1], we define 
the norm of such a lattice L, with respect to H (or equivalently A), as the 
two-sided Oo-ideal i^AiL) generated by the set {H{u,v) : u,v £ L}, and its 
dual by 

L* := {n G L>2 : H{u, L) C Od} ■ 

We see that va{L)"^L C since 

H{au,v) = aH{u,v) € va{L)~^va{L) = Od Va E va{L)~^, u,v E L. 

We say that L is modular with respect to H (or A) if there is an equality, 
and that L is integral if the restriction of to L x L is C'£)-valued. We 
recall that a lattice is said to be maximal if its is maximal with respect to 
the usual inclusion among all lattices of the same norm. We see that if L is 
integral and maximal then it is self-dual, hence modular, since VAiL) = Od- 
From this and the classification provided by [21', Propositions 4.5 and 4.6], 
it follows that all the integral maximal lattices belong to the same genus, 
which we call the principal genus of H. 

Let L and M be two C'£)-lattices in D^. We recall the definition of the 
index [L : M] of M in L. First, assume that M <^ L. Then, the quotient 
L/M is a finite C^-module which admits a composition series 

L/M = Lo D Li D ■ ■ ■ D Lr D Lr+i = 0, 

in which Li/Li^i is isomorphic to Oo/T^i for some maximal left Oc-ideal. 
By the Jordan-Holder Theorem, the set rrij is uniquely determined up to 
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isomorphism by L/M. We define the index of M inside L by 

r 

[L:M] := J]nr(mi). 

1=0 

In general, one defines the index of M in L as 

[L : M] := [L:Lr\M][M : M n L]~^ = [L : P][M : 

where we choose any Oi^-lattice P G LCi M. For the proof of the foUowing 
lemma we refer to Coulangeon [H Lemma 2.2.1], where one can easily see 
that the base field Q plays no special role. 

Lemma 1. Let L = aici © 0262 be an Oo-lattice, where ai and 02 are 
fractional Od -ideals. The index of L in is a square and we have 

[L* : L] = (detD(i^(ei,e,))nr(ai02))^ 

We define the discriminant of L to be the square root of [L* : L\. We 
then obtain the following lemma. 

Lemma 2. Let L and M he lattices in D^. Then, we have = bL[L : M]. 

From Lemmas [Hand EJ it follows immediately that one can equivalently 
define [L : M] as the Oi^-ideal generated by the set {det£)(aa*) : Lot C M}. 
And we see that b^ = i^^(L)^ when L is modular. 

Lemma 3. Let L and M he two OD-lattices in which belong to the same 
genus. Then VA{Mf = va{LY[L ■ M]. 

Proof. It is enough to prove this equality locally. So, let p C Oi? be a 
prime. By definition of genera, we have Mp = L^gp for some gp G GL2(-Dp) 
such that gpAcjp = va{9p)A with VAidp) G -^p^- This implies that VA{Mp) = 
h'A{Lp)vA{gp) and deioigpgl)'^ = det(stp5*) = h'AigpY- Hence, we have 

[Lp : Mp] = (deti5(<7p5p)) = (z^a(5p))' = VA{MpfvA{Lp)-^ . 

□ 

Proposition 4. Let L and M he maximal OD-lo-ttices in such that L 
is modular. Then M belongs to the genus of L if and only if M is modular 
and the index [L : M] is a square Op -ideal. 

Proof. An easy local calculation shows that modularity is a genus property. 
But for a modular lattice M, we have i^AiMf = va{LY[L : M] by LemmasE] 
and[3j Thus we can write ua^M) = UA{L)b, where b is a two-sided O^-ideal. 
The classification of genera of maximal lattices given by [2H Propositions 
4.5 and 4.6] then implies that M belongs to the genus of L if and only if b 
is an C'iT'-ideal; i.e., if and only if [L : M] is a square ideal. □ 

Corollary 5. Let L he a maximal OD-lattice in . Then L belongs to the 
principal genus of H if and only if L is modular and is a square. 
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Let L be a left 0£)-lattice. Since Od is maximal, it is hereditary; hence by 
Reiner [I9l Theorem 27.8], there exists a left Oij-ideal o and a € GL2{D) 
such that L = {On © a~^)a. Furthermore, by the strong approximation 
theorem (or [22^ Theorem 6.14]), we can choose o to be two-sided since the 
class number of (the maximal order) M2(Od) is the narrow class number of 
F. The following generalizes Hashimoto-Ibukiyama |13l Proposition 22]. 

Proposition 6. Let L = (C'D©a~^)a be a maximal On-lcLttice in , where 
a is a two-sided ideal and a G GL2(-D). Then L belongs to the principal genus 
of H if and only if aAa*' = mB, where m £ I'AiL) n and 

B e GU{D) n(^^ ^ . 

^ ' \ a nr(a) 

is a totally positive hermitian matrix such that the ideal detDiB)ni{a)~^ is 
a square. 

Proof. We first observe that, by Lemma [U we have 

(detz)(aa*)detz)(A)) 
nr(a) 

Suppose that L belongs to the principal genus of H. By Corollary [5l this 
means that L is modular and b^^ = i^yi(L)^ is a square ideal or, equivalently, 
that = i'a{L)^^L and yA{L) is an C'jT'-ideal. This implies that 

u{aaAa^)v^ = {au){aAa^)v^ £ OoyaG VA{Ly^ f] F, u,v £ Od ® 

Therefore, we must have aAa^ = mB for some m £ i^AiL) H F^^ and 

'Of a 



B G GL2(-D) n 



nr a 



We see that B is totally positive; and the discriminant formula and the 
relation = va{LY show that det£)(i?)nr(a)~^ is a square. 

Conversely, if there exists a pair {m,B) as above, then b^ is a square. 
Thus we only need to show that L is modular. To this end, we write L = 
OdGi © o~^e2 and L* = Ooef © ae* where the basis {61,62} is given by 
the rows of a. In that basis, the hermitian form is given by the matrix 
aAa^ = mB where 

\r t J ^ ' \ a nr(a) 

Using the facts that m G i^AiL) H F^^ and det£i(-B) G nr(a), we see that 



mdetoiB) \-r s J nr(a) \ a Op 

ML)-' 



Of 

0^-*^ nr(a)^^ 



By combining this with the equality {ef,efY = (aAa^) ^(61,62)*, we get 
the inclusion L"^ C ija{L)~^ L; hence L is modular. □ 
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Remark 7. Proposition [6] provides a simple criterion to search for lattices 
which belong to the principal genus of H; and in practice we can assume 
that m = 1. Indeed, for every m £ F^'^, p2l Proposition 6.13] shows that 
there exists 6 G GL2{D) such that 6A6^ = mA, thus we can replace a by 
a6~^ if necessary. Furthermore, it is not hard to see that we can choose B 
such that detD{B) is a nr(a)-unit. 

From now on, we will denote the lattice L = {On © a^^)a by La^a and its 
discriminant by ba,o. 



1.2. Inner forms of GSp4 and integral models. Let A he a totally 
positive hermitian matrix. We obtain a reductive group G^/Q by putting 

G^iR) := GV^iF 0q R) = {j£ M2(D ®q R) " ""^^^^^ 



Ml) ^{F'^qRV 

for any Q-algebra R. This determines a morphism ua '■ G^ Gm of alge- 
braic groups called the similitude factor. The group G'^/Q is an inner form 
of Resp'/(Q(GSp4) and, by combining [21, Proposition 2.1] and [lOl Lemma 
4], we can obtain an explicit isomorphism G^{E) ~ GS'p^{Ey . 

Let L be a maximal integral O/j-lattice in D^. We obtain an integral 
structure on G^ over Z by putting 

7^7* = '^a{7)A 
UAi-f) G {Of ®i.RY 



(i?) := |7 E Endo^552R(L ®z R) 
= G\]^{L®i R), 



for any Z-algebra R. We call the principal genus of G^ , the one containing all 
the integral models obtained from integral maximal lattices in the principal 
genus of H. We observe that, for every prime p, 

GiiZp) = llG\J^iL Of,), 
p\p 

where G\J^{L ©Fp) is hyper special if and only if p ^ S. Therefore, G^ 
is an integral model in the sense of Gross [T2] if and only if S = 0. This is 
equivalent to saying that [F : Q] is even and D/F is the unique quaternion 
algebra ramified at | oo only, in which case there is only one genus. 



2. Algebraic modular forms and Hecke action 

From now on, we let A = I2 and we simply denote the group G^ and its 
integral model associated to L = Ojj by G and G, respectively. We recall 
the theory of algebraic modular forms on G in the sense of Gross |11) . We 
then introduce the notion of algebraic modular forms on integral models of 
G and explain how the two relate. 
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2.1. The Hecke module of algebraic modular forms. For any finite 
prime p ^ S, we choose an isomorphism GU2(ODp) — GSp4(C'Fp), which is 
compatible with the splitting isomorphism O^p = ^2{Of^ ) we fixed earlier, 
and we consider the maximal compact open subgroup 

G{Z) = H GV^iOn,) X n GSp4(OfJ. 

pes p^s 

We let U = Hp Up C G{Z) be a compact open such that Up is maximal for 
each prime p G S. 

Let p : G ^ GL(y) be an irreducible algebraic representation. For any 
subfield Q C C C, we denote the realization of (p, V) on E in the same 
way. We recall that, when E splits D, then this realization is of the form 
P = 'S)v\oo P'" V = <S)v\ao ^v, where each 

: GSp4(S) ^ GL(K) 

is an irreducible representation of highest weight (ai,a2; b) G Z^, with ai > 
02 > 3 and b = ai + a2 mod 2. 

The space of algebraic modular forms on G(Q) of weight V and level U 
is defined by 

MpiU) := {/ : Gm/U ^ V\ f\pj = /, for ah 7 G G(Q)} , 

where f\p-f{x) = f{jx)p{-f), for all x E G(Q) and 7 G G(Q). The Hecke 
algebra Tp{U) acting on this space is the Z-algebra generated by the oper- 
ators defined as follows. For any u G G(Q), write the finite disjoint union 
UuU = \JiUiU and put 

f\p[UuU]{x)=Y,fi^U^)■ 

i 

For a prime p ^ S, the local Hecke algebra at p is generated by the two 
operators Ti(p) and T2(p) corresponding to the double cosets of the diagonal 
matrices diag(l, 1, -ojp, vjp) and diag(l, -ojp, zup, tUp) respectively, where Wp 
is a uniformizer at p. Hence the Hecke algebra Tp(U) acting on Mp(U) is 
the Z-algebra generated by the operators ri(p) and T2(p) for all p ^ S. 

When p = triv is the trivial representation, we let Ip{U) be the subspace 
of Mp{U) given by 

Ip{U) := {/ : G(Q)\G(Q)/f/ ^ FX+\(A^)x/i.(C/) c} . 

Then, we define 

r Mp{U), p/triv 
Sp{U) := <^ 

[ Mp{U)/Ip{U), p = triv. 

The action of Hecke preserves Ip{U), and thus induces an action of Tp{U) on 
Sp{U). We refer to lhara, Hashimoto- Ibukiyama, and Ibukiyama [TBlflSlfTS] 
for the conjectural analogue of the Jacquet-Langlands correspondence for 
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GSp4 (over F = Q) which relates the Hecke module Sp{U) to spaces of 
Hilbert-Siegel modular forms, and in particular to Sorensen |23j who shows 
that this conjecture is true when the degree of F is even. 

2.2. Integral models of G and the Hecke module Mp{U). Let C1(G) 
be a complete set of representatives of the conjugacy classes in the genus 
of G as described in p^. There is a canonical bijection between C1(G) and 
the double coset G{Q)\G{Q)/G{Z), which also parametrizes isomorphism 
classes of left Ozj-lattices in the genus of Ojj. By abuse of notation, we let 
(a, a) be the integral model attached to the lattice -La,a- Let Sf/ := (S^, T,fj) 
be a pair of sets consisting of primes such that n = and each is 
disjoint from T, and the set of primes at which Up is not maximal; further 
assume that generates the narrow class group Cl^(-F). By the strong 
approximation theorem we can choose the representatives (a, a) £ C1(G) 
to be supported in T^u; i.e., such that nr(o) is supported in and a is 
S^-integral. (In practice, one can choose the set to be a singleton, but 
it can't be empty unless the class number of G is equal to := ^CV' (F) .) 

We define J^cjU) '■= G{'Z)/U (this is a finite product of finite sets). Let 
(a, a) S C1(G) be a representative; the group of automorphisms of the cor- 
responding lattice La,a is given by 

Ta,a := StabG(Q)(La,a) = S^l & G(Q) D a'' (^^ a : G 

and the quotient Ta,a = ^a,a/C)p is finite. We define 

Mp""(f/) := {/ : ^FciU) ^ V : /|p7 = / V7 G r„,„} . 
Let (b,/5) G C1(G) be another representative, and recall that 

Homo^(L„,„, Lfa,^) := {7 G M2{D)\ L,,„7 C Lf,,^} = (^^^ /3. 

Now, let p ^ S be a prime ideal in Oi;', and choose a' G a^^ and 6Gb 
such that Op"*^ = (a'p) and bp = (6p); (this is equivalent to saying that 
fp(nr(a')) = — Wp(nr(a)) and Vp(nr(6)) = Wp(nr(b))). Then, for every element 
7 G Homc)^(La^a, Lb,/?), the matrix diag(l, a')a^P~^diag{l, b) is p-integral. 
Thus, we get a well-defined map 

^(a,a),(b,/3) : Homci^(L„,„, Lt,,/?) ^ M2{Od,) 

'I 0\ 



^ ^ (0 a°p)"P>^l 



P \0 bp 

Let Fp = Of/P be the residue field at p and Od^, ^ 1^2(0^^) ^ M2(Fp) a 
reduction map (obtained from a fixed splitting map). Then, by composition 
with i?(a,o),(t),/3)5 we get a well-defined map 

Homo^(La,„, Lb,/3) ^ M4(Fp) 

71 — > 7 
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Next, for i = 1 and 2, we define the sets 
ef^(p;(a,a),(b,/?)) := 

Ta,a\ e G(Q) nHomo^(La,a, Lj,,/? 

where we let Ta,a act on the left. By unravelling definitions, we get 



[Lb,f3 ■■ La,a7] = P\ 
rankp (7) = 3 - i 



G^(p;(a,a),(b,/3)) := 



ra,a\ <^ 7 G 



n a" 



-1 (Od b-i 



ah' 



1 /5 I^(7)fca,Qbf,_^ 



e2^^(p;(o,a),(b,/3)) := 



n a 



1 



a ab 



1/3 



I^(7)fca,abt,,J = 
rankFp (7) = 1 



(It is not hard to see that these sets do not depend on the choice of a' G a ^ , 
6Gb and the splitting map at p.) We define the linear map 



r,(p;(a,a),(b,/5)) : M;'-'{U) 

f 



/Ipn, i = l,2. 

«eef^(p;{a,a),(b,/3)) 

The following proposition is proved in the same way as [6^ Theorem 2]. 
Proposition 8. There is an isomorphism of Hecke modules 

(a,a)eCl{G) 

where the action of the Hecke operator Ti(p) on the right is given by the 
collection of linear maps Ti{p; {a, a) , {b , P)) , (0, a),(b,/3) G C1(G). 

Proposition [8] also gives a decomposition ^^(C/) ~ 

(a,a)eCl{G) 



3. The algorithm 

We now explain how to implement the discussion of Section [2] into a 
concrete algorithm. As with our presentation, the algorithm we present 
here parallels the one in [5] and [6]. 

3.1. Representing hermitian quadratic forms. In this subsection, we 
discuss the problem of representing hermitian quadratic forms by hermitian 
quadratic forms. This problem arises naturally when testing for equivalence 
in C1(G) or during the computation of the Hecke action. We remark that 
analogous questions appear when dealing with Fourier expansions of Siegel 
modular forms. 
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Let j4 be a totally positive hermitian matrix and Q the associated qua- 
dratic form. The matrix A also gives a vector- valued Z)-hermitian form 

M2{D) X Map) ^ M2{D) 
(71,72) ^ 71^72- 

Let Oi, bi, Ci and Oi be left O/^-ideals in D so that A := {^^ | is an 

O/j-lattice in M2(D). Given a hermitian matrix 77 € M2(-D), we want to 
compute the set 



0vl(r/) := <^ 7 G A : = ?? 



s r 
r t 



Let 7=1 I G M2(-D) be given in rows. Then, it is easy to see that 



7A7* 



s r 
r t 



{ Q{u) = s (1) 
vAu^ = r (2) 
Q{v)=t (3) 



Thus computing 0a(??) is the same as representing s,t G F by Q under the 
constraint (2). This implies that 0A(f?) is finite since ^ > 0. By Shimura [21^ 
Proposition 2.1], there is a diagonal matrix r/' G M2(-F) and P £ GL2(-D) 
such that rj = Prj'P^. By making the change of variables 7 = P7' and letting 
A = PA', we see that this is again the same as computing the set 

Qa{v') := {7' € A' : y^f* =V'= (q °) } , 

where the constraint (2) is now an orthogonality condition. Let A'[l] and 
A' [2] be the rows of A', which are Oj^-lattices in D^. To compute 0a(?7')) 
we first determined which of the row lattices (A'[l], Q) or (A' [2], Q) is better 
to work with by comparing their determinants. Without loss of general- 
ity, assume this would be A'[l]. We then apply the algorithm in Dembele- 
Donnelly [3 Subsection 2.2] to efficiently compute the set of elements in 
(A'[1],Q) which represent s' . 

Now, let u = {a,b) G A'[l] be a solution to Q{u) = s' . We replace u 
into (2) to get a homogenous linear system in terms of the components of 

V = (c, d) G A' [2] . Solving this linear system in terms of c or d and then 
replacing into (3), we must now represent t' by a quadratic form in 
one variable over an Oo-^desd. Prom each solution of Qu{w) = t' , we get 

V £ A'[2] which, combined with u, gives a matrix 7' G QAiv')- (The trick of 
using the orthogonality condition (2) to cut down the number of variables 
in (3) is substantially more efficient than the naive approach of solving (1) 
and (3) separately and then checking (2) for compatibility.) 
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3.2. The representatives (a, a) G C1(G). We explain how to compute 
representatives for the equivalence classes in the genus of G. In order to do 
so we need the following two lemmas both adapted from |13j . 

Lemma 9. Let L^^a o-nd La^^i he in the genus of G. Then, L^^a cLnd La,a' 
are equivalent if and only if there exist n G F^~^ and (3 £ GL(C'/) © a^^) 
such that a'a'^ = nf3{aa^)[3^ . 

Proof. This easily follows from Proposition [6l □ 
Let a be a two-sided ideal. Then, motivated by Lemma [9l we say that 



there exist n E and 13 £ GL{0 d ® oT^) such that 7' = n(5-f(5^. 

Remark 10. We test whether two lattices L^^a and L^^p in the genus of 
G_ are equivalent by first reducing to the situation in Lemma [H By the 
observation preceding Proposition [6l this amounts to testing whether a and 
b are in the same narrow class. In fact, Proposition [6] and Lemma [9] imply 
that the equivalence class of the lattice L^^a is completely determined by 
the narrow class of a, and the a-equivalence class of 7 := aa* mod F^'^ . 
(See Remark [71 and [211 Lemma 4.4] by which we can always find a matrix 
a G Gh2{D) such that aa* = 7.) 

Lemma 11. Let a be a two-sided Od -ideal and s G O^, and let 



Then, the a-equivalence class of ^ depends on the class r mod sa only. 
To compute the set C1(G), we use the following algorithm. 

Algorithm 12. Given a finite set of primes T, and an integral ideal n sup- 
ported outside of S, this returns a set C1(G) of representatives (a, a) sup- 
ported outside of S and n. 

(1) Initialize Rq = and Mq = 0, and compute the mass Mass(G). 

(2) Find a set S' of prime ideals that generate Cl~''(-F) and such that S' 
is disjoint from Sujp : p | n}. Choose a prime q ^ S'UllUlp : p | n}. 

(3) Compute a set Sd of integral two-sided O^i-ideals supported in S' 
and such that, for each b G Cl^(F), there is a unique a G Sd such 
that [nr(a)] = b^. 

(4) For each a G Sd, find a set of representatives Ua for the quotient of 



[ nr(a) J 

modulo squares. (This is equivalent to a Selmer group computation.) 

(5) Choose s G Op (ordered by norm). 

(6) Choose a G Sd and compute the set a/sa. 



two hermitian matrices 7, 7' G GL2(-D) H 



Of a 
a nr(o) 



) 



are a- equivalent if 
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(7) For each {u,r) £ Ua x such that nr(r) + u G nr(a)s, put 



is aheady represented in Rq- If not, solve the equation aa* = 7 for 
a matrix a, which is q-integral, and compute the finite group Fa^a- 
Append (a, a) to Rq, set Mq = Mq + ^f— • 



(8) If Mg = Mass(G), then C1(G) = hence stop. Else, go to Step 
6. If Step 6 fails, go to Step 5. 

Remark 13. We test the equivalence of two O/j-lattices L and L' by using 
Lemma [9] or by comparing their theta series. The latter is often enough to 
quickly distinguish most isomorphism classes and we only use Lemma [9] for 
the remaining cases. 

Let -L be a modular lattice so that ^{L) is an Oir-ideal, and let c be a 
representative of the narrow class of u{L)^^ . We choose S such that 
u{L)c = (fj,). For every element (x,y) G L, we define its scaled norm with 
respect to L and c by 

A* 

For a given representative c, this is well-defined up to a totally positive unit. 
For every integral ideal m C Op, we define 



where F^, := AutGr(Q)(L). Then -^^(m) is non-empty if and only if ^{L)^^ 
and m are in the same narrow class. In that case, let it be a totally positive 
generator of u{L)m. Then the set of elements (x, y) £ L such that i'{x, y) = u 
is clearly finite. Thus the quotient Rii^) is finite. Let rL{m) := #i?L(m) be 
its cardinality. For each t £ Sj^ , where Sj is the Poincare upper-half plane, 
we define 



where the second summation is restricted to all m C O^? that belong to the 
narrow class of c, and for each such m we let i^c(m) be a totally positive 
generator of mc~^. A similar argument as in Eichler [HI &: 4, Theorem 1] 
yields that is a Hilbert modular theta series of parallel weight 2 and 
whose level depends on c and the primes in S. 

3.3. The stabilizers F^.q and the sets Gj ^ {p; {a, a), In this sub- 

section, we describe how the stabilizers and the prelimimary data required 
for the Hecke action are computed. To this end, let (a, a),(b,/3) G C1(G) 



t := (nr(r) + u)/s. Check if the a-equivalence class of 






{x,y)GL/rL 



mCOp 
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and u ^ F totally positive. We can compute the finite set 

S{u; (a, a), (b, /?)) := {7 G G(Q)| La,a-/ ^ L^^p and 1/(7) = u] 

by making use of the algorithm in Subsection 13.11 where we first write out: 
_i iOd i>'^ \ o_ ^ai t'l 



P\ 77* = 



a 



ci 



To obtain T^^a (or equivalently Ta,a), we simply set (a, a) 

ives for 

]J S{u;{a,a),{a,a)). 



let Up be a set of representatives for Op~^/ (Op) . Then, we get 



[b,/3) and 



For the computation of the representatives of the Hecke double cosets, we 
first observe that, for i = 1,2, the set 0j ^ {p; (0, a), (b,/3)) is non-empty if 
and only if the ideals bajob^^ and are in the same narrow class. Letting 



u 



be a totally positive generator of H./J^a «P*' then get 

ef ^ (p; (0, a), (b, /?)) = r,,„\5(n; (a, a), (b, /?)). 

Remark 14. For each prime p C Op, the cardinalities of these sets satisfy 
the following identity (see [71 Remark 3]): for i = 1,2 and (b,/3) E C1(G), 

^ #ef^(p; (a, a), (b,/3)) = Np-i(Np + l)(Np2 + 1) = degT,(p). 

(a,a)eCl(G) 



3.4. The description of the flag spaces Tg{U). In this subsection, we 
describe the flag J-g{U). The level structure U will be one of the following 
types: Siegel, Klingen or Borel parahoric (at primes p ^ S). Let n C Op be 
an ideal such that p f n for p G S. 

• The Siegel parahoric 



U :-- 



7 G GiZ) : 7 ^ 



* 
* 



mod n 



7 



* * 

\o 

We recall the following description from [5] . Consider the free rank 4 (Op/n)- 
module E := (Op/n) endowed with the symplectic pairing (, ) given by 
the matrix 

where I2 is the identity matrix in M2{Op/n). We fix the canonical sym- 
plectic basis ei,e2,/i,/2 such that (cj, ej) = {fi, fj) = and (cj, /,) = 5ij 
(Kronecker delta). Let M be a rank 2 (C'ir/n)-submodule of E. We say that 
M is isotropic if {u, v) = for all u, v & M. Let .^^(n) be the set of all rank 
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2 (Oi?/n)-submodules of E which are isotropic direct factors. We recah that 
G(Z) acts transitively on J^ci^), and that the stabiUzer of the submodule 
generated by ei and 62 is the Siegel parahoric U. In other words, J^g{^) is 
the set of (C'ir/n)-rational points of the Lagrange scheme. 

By using Pliicker's coordinates, we can identify J^ci''^) with the set of 
(Oi^/nj-rational points of the Klein quadric as a closed subspace of P^(C'i7'/n). 
We then represent each element in J-g{t^) by choosing a point x = (oq : • • • : 
as) = [« A G P^(C'ir/n) such that the submodule M generated by u and 
u is a rank 2 isotropic direct factor of E. We observe that the module M 
is a direct factor of E if and only if the ideal (ao, ... ,05) generated by its 
coordinates is the unit ideal. 

The cardinality of .^^(n) is given by [5, Lemma 3.2]: 



#J^G(n) = N(n) 



n 

p|n 



The Klingen parahoric 



U :-- 



7 e G{Z) : 7 ^ 



1 + 



* 



1 



N(p) 



1 + 



N(p)2 



mod n 



> . 



Consider the module E as above and let .^^(n) be the set of all rank 1 
(C'ir/n)-sub modules which are direct factors of E. The group G(Z) acts 
transitively on J-ci^) and the stabilizer of the line L = (ei) is the Klingen 
parahoric U. This means that we can identify Tcin) with the projective 
space P^(C'i?/n). Thus, its cardinality is given by 

#J^G{n) = N(n)3 n f 1 + Jr^] (l + ^ 
p|n 



N(p) 



N(p)2 



The standard Borel parahoric 

/* 

U := (-fG : 7 = 



* 





* 





*/ 



mod n 



Again, we let E be the module as above. We let Tg{^) be the set of all pairs 
{L,M) such that L C M are rank 1 and 2 (Op'/n)-submodules of E which 
are isotropic direct factors. The group G{Z) acts transitively on ^g(i^) and 
the stabilizer of the flag ((ei), (ei, 62)) is the standard Borel parahoric U. 

Lemma 15. The cardinality of J-G{n) is given by 

#^e(„) = N(„)'n(i + j4)'(i + i4 
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Proof. Let M be an (Op'/n)-rational point on the Lagrange scheme and 
e'l, 62 a basis of M. A hne L = {aie'i + 0262), with 01,02 S Op/n, is a 
direct factor of M if and only if (oi, 02) determines a point on P^^Op/n). 
Combining this with the formula for the cardinality for the Lagrange scheme, 
we get the lemma by recalling that 



Remark 16. It is useful to observe that the proof of Lemma [l5] gives an 
effective way of listing elements of the flag space of the Borel parahoric. 
Indeed, by using the description of the Lagrange scheme, we can efficiently 
list its elements. Then for each such element, we list all the lines that are 
direct factors by simply running through P'^i^Op/n). 

Remark 17. Lansky and Pollack [T7] already make use of the flag space 
Tg{U) in order to find coset representatives for the double coset space 
G{Q)\G{Q)/U. But, they then return to the adelic setting for the com- 
putation of the Brandt matrices. One of the advantages of our approach is 
that, using Proposition [H we define the Hecke action on Tg{U) directly. 

3.5. A sketch of the algorithm. Here, we give a sketch of the implemen- 
tation of the algorithm. 

Algorithm 18 (Precomputation). The input is a field F as above, a finite 
set of finite places S such that + [F : Q] is even, and an integral ideal n 
not supported in S. 

(1) Find a quaternion algebra D/F ramified at precisely the infinite 
places and the places in S, and compute a maximal order Od of D. 

(2) Find a set E' of prime ideals not dividing n that generate Cl^(-F) 
and such that S n S' = 0. Choose a prime q ^ S U S' U {p : p | n}. 

(3) Compute a complete set C1(G) of representatives (0, a) for the classes 
in the genus of G which are supported in (E', {q}) (see Algorithm [12]). 

(4) For each representative (a, a) G C1(G), compute the group Ta^a. 

Algorithm 19 (Main algorithm). The input consists of F, the level n (with 
the level structure specified: Borel, Klingen or Siegel) and the precomputed 
data, together with a weight representation (p, V). The output consists of a 
set of primes p (ordered by norm), the Hecke matrices Ti(p) and T2(p) and 
the Hecke constituents. 

(1) Compute splitting isomorphisms 0X12(00,) — G'Si\>^{Of^), for each 
prime q | n. 

(2) For each (a, a) G C1(G), compute Mp'''°(n) as a module of coinvari- 




□ 



ants 



Mp»-(n) = E[TGixx)\ ®V/{x- 7X, 7 G r„,,). 
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(3) Combine the results of step (2), forming the direct sum 

M,(n)= Mp^-(n). 

(a,Q!)6Cl(G) 

(4) Starting with the primes of smallest norm in Op, compute the Hecke 
operators ri(p) and r2(p) for all the primes p, in increasing order, 
until you find a common basis of eigenvectors of Mp(n) for the ri(p) 
and T2(p) that completely diagonalize Mp(n). 

Remark 20. In practice, one can improve the efficiency of Step (4) of 
Algorithm 1191 by diagonalizing Mp(n) only using the Hecke operators ri(p) 
and then computing the T2{p) when needed. 



4. Examples of Hilbert-Siegel eigenforms on 

In this section, we give some examples of Hilbert-Siegel eigensystems for 
F = Q(-v/2). (We obtained them with a preliminary version of our algorithm, 
which we implemented in Magma [2j.) 

4.1. Hilbert-Siegel eigensystems of level 1. Let D be the Hamilton 
quaternion algebra over F, i.e., the quaternion algebra over F determined 
by the relations = — 1, = — 1 and k = ij = —ji. Then disc(L') = (1) 
and we choose the maximal order Od = ^[V^][ei, ^2, 63, 64], where 

l+i l+j 
ei = —^■< ^2 = ^3 = 6162, 64 = 6261. 

We recall that the narrow class number of F and the class number of the 
quaternion algebra D are both equal to 1. Thus every representative in the 
principal genus of G is of the form ((1), a), where aa* = 7 G GL2(0_d). We 
choose = (0, {(\/2)})- We determine that the matrices 

_ /I 0\ _ / 2 (2W2)+(2-V2)l\ 

~ 1^0 1 y ' ^2 ~ 1^ 2+x/2+(-2+72)» 2 j 

correspond to distinct equivalence classes by computing their stabilizers for 
which we obtain Card(ri) = 4608 and Card(r2) = 3840, where we simply 
let Tj = r(i) Q. . By using the mass formula in W. T. Gan, J. Hanke and 
J.-K. Yu Proposition 9.3], we verify that 

1 1 1 11 1 1 

Mass(G) = _Cq(^^(-1)Cq(^)(-3) = ^ — — = ^ + ^ 

1 1 

+ 



Card(ri) Card(r2) 

Therefore, the class number of the unique genus of G is h = 2. 

This means that the dimension of the space of holomorphic Hilbert-Siegel 
modular forms of level 1 and parallel weight 3 is 2. From this, we determine 
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that there is only one Hilbert-Siegel cuspidal eigenform / of level 1 and par- 
allel weight 3. The Brandt matrices 13i{p) and ^2(p) of the Hecke operators 
acting on M^^l) for the first four prime ideals p are given below. 



N(p) 



{V2) 



7 (3 + ^2) 



7 {3-V2) 



(3) 



Blip) 



^2(P) 



9 6 
5 10 



208 192 

160 240 

208 192 

160 240 



12 18 
15 15 



1264 1536 

1280 1520 

1264 1536 

1280 1520 



436 384\ /3540 3840 
320 500y 1^3200 4180 



The first few Hecke eigenvalues of the eigenform / are listed in Table [TJ 
The eigenvalues of / (or in fact the Brandt matrices) at the primes above 7 
suggest that it is a base change from Q, which we then confirm as follows. 
For a classical or Hilbert newform h of weight 4, let SK{h) be its Saito- 
Kurokawa lift. Letting ap{h) be the Hecke eigenvalue of h at p, the Hecke 
eigenvalues Ai(p) and A2(p) of SK{h) corresponding to the Hecke operators 
Ti(p) and r2(p) are then given by 



Ai(p) 

A2(P) 



ap{h) N(p)^+N(p)2 + N(p) 
a^ih) N(p)^(N(p) + l) + N(p)2 



1. 



Proposition 21. The form f is a base change of the form SK{h), where h 
is a newform which belongs to the unique conjugacy class o/ 54(8, (-^))"'^™, 
the newspace of classical forms of level 8 and character {-p). 

Proof. We compute the space 54(1) of Hilbert cusp forms of level 1 and pa- 
rallel weight 4 over F using the Hilbert Modular Forms Package in Magma [2], 
and obtain that it has dimension 1. We then identify the eigenform / as the 
Saito-Kurokawa lift of the newform g in 54(1) by direct calculations. By 
Saito j20| Theorem 4.5], the form g itself is a base change of a newform h in 
54(8, (7))"™. The hft SK{h) is a classical Siegel eigenform of level 8 and 
character (-p). By the funtoriality of Saito-Kurokawa lifts |3J, we get 



/ = SK{g) = SK{BCp/Q{h)) = BCF,Q{SK{h)). 



□ 



Remark 22. Proposition [211 combined with Sorensen [23, Theorem B], 
implies that there is no holomorphic Hilbert-Siegel eigenform of level 1 and 
parallel weight 3 over Q(\/2) which is stable. 
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By lifting of Galois representations attached to /i, there is a family of 
^-adic Galois representations ipsK{h),e) attached to SK{h): 

PsK(h),i ■■ Gal(Q/Q) ^ GSp4(Q,). 

Similarly, there is a family (p//) of ^-adic representations of Gal(Q/F) at- 
tached to /, which is everywhere unramified since / has level 1. In fact, by 
Proposition EH we simply have pf^i = PsK{h),e\Gs.m/F)- 

It would be interesting to know whether there exists a (Calabi-Yau) three- 
fold X/q{V2), with Hodge numbers = /i^'^ = /i^'^ = /i^'S = such that 
for every prime £ the Galois representation 

Px/ : Gal(Q/Q(^/2)) ^ GL(i/3 (X, Q,)) 

is isomorphic to pf^e. The existence of such a threefold would provide an 
analogue in higher rank of the Tate (or Shimura) elliptic curve i?/Q(\/29) : 
y'^ + xy + e^y = with discriminant A^; = —e"^^, where e = ^^^^ is 
the fundamental unit, in the sense that Q(\/2) would be the totally real 
field of smallest discriminant for which there is a higher rank (automorphic) 
motive over Q whose base change to Q(\/2) becomes everywhere unramified. 
(Indeed, the results in [5J show that no such threefold exists over Q(\/5)-) 

4.2. Hilbert-Siegel eigensystems of prime level. For all prime levels 
r, with odd norm N(r) < 31 (and up to Galois conjugation), we compute 
the space Sp{x) using the same precomputed data as in 14. H for the triv- 
ial representation p and the Siegel parahoric. By the Jacquet-Langlands 
correspondence for GSp4 for fields of even degree [23l Theorem B], we ob- 
tain Hecke eigensystems of Hilbert-Siegel cuspidal forms of parallel weight 
3. (Note that [23^ Theorem B] is only true for eigensystems which come 
from stable automorphic forms; but for forms that are lifts, one can use the 
usual Jacquet-Langlands correspondence for GL2.) In Tables [T] and [21 we 
list all those eigensystems that are defined over Q or a quadratic field. Here 
are the conventions we use in the tables. 

(1) For a quadratic field K of discriminant d, we let uod be a generator 
of the ring of integers Ok of K. 

(2) The first row contains the level r given in the format (N(r), a), where 
a € Op is a generator oft, and the dimensions of the relevant spaces. 

(3) The second and third rows list the computed Hecke operators. 

(4) The Hecke eigenvalues systems are given by row. 



r = (1, 1), dimM3(r) = 2, dimS'3(r) = 1 


N(p) 2 7 


7 




9 




Ti{2 + co8) T2{2 + U8) Ti{3 + i08) 12(3 + ^8) Ti(3 


-tJs) 72(3 


- ^s) 


Ti(3) 


^2(3) 


/i 4 -3 48 -16 


48 


-16 


116 


340 



r=(7, 3 + a;8), dimM3(r) = 


6, dimS'3(r) = 


= 5 














N(p) 2 


7 








7 






9 




ri(3 + tc;8) 


r2(3 + a;8) 


ri(3 


- t^s) 


T2(3 


- W8) 


Ti(3) 


7^2(3) 


h 10 15 


-7 







60 




80 


80 


-20 


/2 -4 1 


7 







32 




80 


-60 


148 




r=(9, 3), dim Ma (r) = 12, 


dim 53 (r) = 11 














N(p) 2 


7 








7 






9 


ri(2 + a;8) T2(2 + a;8) 


ri(3 + a;8) 


72(3 + a;8) 


ri(3 


-UJs) 


72(3 


-<^8) 


ri(3) 


r2(3) 


/i -6 7 


-22 


64 




-22 




64 


-9 





/2 4 5 


-16 


16 




-16 




16 


9 







r = (17, 5 + 2u;j^), dimi\/5(c) 


= 23, dim Sy, 


(r) = 22 














N(p) 2 


7 








7 






9 


ri(2 + W8) 12(2 + u;8) 


ri(3 + w8) 


72(3 + 1^8) 


Ti(3 


- Ws) 


72(3 


-U8) 


Ti(3) 


7^2 (3) 


/i -2 3 


6 


32 




-36 




80 


8 


28 


/2 4 - a;4o -3 - 3a;4o 


38 - Aujio - 


96 — 32a;4o 


58 + 6a;4o 


64 + 48a;4o 


66 + 4a;4o 


-160 + 40a;4o 



Table 1. Hilbert-Siegel eigenforms of parallel weight 3 over Q(\/2) 



r= (23, 5 + a;8), 


dimM3(r) = 


32, dim53(r)=31 










N(p) 


2 


7 


7 




9 






T2(2 + W8) 


Ti(3 + W8) 12(3 + (jg) 


ri(3-u;8) T2(3 


-Us) 


Ti(3) 


^2(3) 


fi -8 


10 


-46 119 


-54 


145 


-54 


153 


f2 10 


15 


72 176 


56 


48 


124 


420 



r=(25, 5), dimM3(r) =48, 


dim 53 (r) = 47 














N(p) 2 


7 




7 






9 




ri(2 + a;8) T2i2 + ujs) 


Ti(3 + (J8) T2(3 + W8) 


ri(3 - ujs) 


r2(3 - LOs) 




ri(3) 




^2(3) 


h 2 -9 


62 96 


62 


96 




40 




-420 


/2 -4 - 2u;i2 5 + 2a;i2 


— 15 + 7ix'i2 56 — 8u;i2 


-15 + 7a;i2 


56 — 86<Ji2 


44- 


18wi2 


132- 


- 52t<;i2 


/3 8 9 


44 - 10a;24 -48 - 80u;24 


44 + 10a;24 


-48 + 80u;24 




76 




-60 




r = (31, 7 + 3a;8), dimM3(r) 


= 65, dimS3(r)=64 














N(p) 2 


7 




7 






9 




ri(2 + a;8) r2(2 + W8) 


Ti(3 + W8) T2(3 + (J8) 




r2(3 - ujs) 




Ti(3) 




^2(3) 


/i -4 4 


-60 174 


-20 


-2 




-60 




130 


/2 -4 5 


20 16 


-8 


48 




8 




28 


h -6 7 


44 - 26J204 112 - 86^204 


— 4^^204 


48 


-32- 


2^204 


108 


+ 46^204 


/4 2 -9 


72 — 2a;204 1 76 — 16aj204 


56 - 4a;204 


48 - 32CJ204 


76- 


2aj204 


-60- 


20a;204 



Table 2. Hilbert-Siegel eigenforms of parallel weight 3 over Q(a/2) (cont'd) 
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